Quantum-classical correspondence in the hydrogen atom in weak external fields 
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The complex processes leading to the collisional popula- 
tion of ultralong-lived Rydberg states with very high angular 
momentum can be explained surprisingly well using classical 
mechanics. In this article, we explain the reason behind this 
striking agreement between classical theory and experiment 
by showing that the classical and quantum dynamics of Ry- 
dberg electrons in weak, slowly varying external fields agree 
beyond the mandates of Ehrenfest's Theorem. In particular, 
we show that the expectation values of angular momentum 
and Runge-Lenz vectors in hydrogenic eigenstates obey ex- 
actly the same perturbative equations of motion as the time 
averages of the corresponding classical variables. By time 
averaging the quantum dynamics over a Kepler period, we 
extend this special quantum-classical equivalence to Rydberg 
wave packets relatively well localized in energy. Finally, the 
perturbative equations hold well also for external fields be- 
yond the Inglis- Teller limit, and in the case of elliptic states, 
which yield the appropriate quasiclassical initial conditions, 
the matching with classical mechanics is complete. 

32.80.Rm, 32.60. +i,34. 10. +x,03.65.-w - Printed in Phys. Rev. 
A 58, 3896 (1998). 



I. INTRODUCTION 

In the past few years new experimental techniques have 
made possible the study of the dynamics of atoms or 
molecules in which an electron is promoted to a very 
high energy state, where it is only weakly bound to the 
core [Q. These high energy states can be described by 
approximately hydrogenic wave functions with very large 
principal quantum numbers (n > 100) |^,^|. The atoms 
(or molecules) in which a valence electron is promoted 
to such high-n states are generically called "Rydberg" 
atoms, because the energy levels of the excited electron 
are well described by a Rydberg-like formula and 
their highly energetic electron is known as a Rydberg 
electron. In such systems the weakly bound Rydberg 
electron resides mostly at an immense distance from the 
atomic core, a distance so large that if Rydberg atoms 
were solid, they would be just about visible to the naked 
eye. Laboratory-scale external fields, and even weak 
stray electric fields , become then comparable to the 
atomic Coulomb field sensed by the Rydberg electron, so 
that the dynamics of the electron can be probed with ac- 
curacy, and also fundamental dynamical properties such 
as quantum manifestations of chaos [p|-|l4|], can be stud- 



ied experimentally. 

To a very good approximation, the dynamics of Ryd- 
berg electrons is hydrogenic. More precisely, small devia- 
tions from the purely hydrogenic eigenenergies are intro- 
duced by the interaction between the far flung electron 
and the electronic cloud around the atomic or molecular 
core. These deviations are quantified by the quantum de- 
fect Si, which enters in the formula for the energy levels 
as a correction to the principal quantum number n 0, [l5| . 
However, the quantum defect becomes rapidly negligible 
as the angular momentum of the electron increases. In 
fact, more complex atoms are often used as experimen- 
tal substitutes for hydrogen, since it is much easier to 
excite their valence electron to a Rydberg state, and yet 
the field sensed by the Rydberg electron does not differ 
much from a pure Coulomb field. Therefore many recent 
investigations of Rydberg electrons in alkali atoms have 
really probed the dynamics of quasiclassical electrons in 
an essentially coulombic potential. These accurate ex- 
perimental results have led to a renewed theoretical in- 
terest in the hydrogen atom in external fields in the limit 
of large quantum numbers, Jl(| [l6] |l9| which has become 
one of theparadigmatic models for the study of quantum 
chaos flyTLJ] , and of quantum-classical correspondence in 
general. 

Since the degeneracy of a hydrogenic n-manifold grows 
as n 2 , a fully quantum treatment of the dynamics of Ry- 
dberg electrons poses formidable challenges even to the 
most advanced computers. Therefore classical mechan- 
ics is often the only practical way to study such sys- 
tems, under the assumption that for large n's classical 
and quantum predictions should somehow converge. In- 
terestingly, however, recent experimental and theoretical 
work on electronic wave packets in hydrogenic systems 
|p0| [l5| have shown very clearly that the quantum me- 
chanical properties of the Rydberg electron are essential 
to the dynamics of the wave packet, even in the large- n 
regime. For example, the observation of fractional re- 
vivals can be explained only by the quantized spectrum 
of the Hamiltonian |20|-|2^, p5| , |27| |29[ | in spite of the very 
large principal quantum numbers involved, which seems 
to suggest that large quantum numbers are not sufficient 
to ensure the accuracy of a purely classical description of 
Rydberg dynamics. 

On the other hand, classical mechanics yields surpris- 
ingly accurate results for the problem of the hydrogen 
atom in weak, slowly varying external fields; i.e., when 
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the classical electron still moves, to a good approxima- 
tion, along a Kepler ellipse, and the semimajor axis of 
the ellipse (or, equivalently, the Kepler energy of the elec- 
tron) remains unchanged. In particular, recent classical, 
perturbative calculations [l6] have succeeded in ex- 
plaining several diverse experimental results with aston- 
ishing accuracy, ranging, for example, from the ultralong 
lifetimes of the molecular high-n states employed in zcro- 
electron-kinetic-energy (ZEKE) spectroscopy §-§f§7|,|4| 
|52|-|58| , to the intrashell transitions induced in alkali Ryd- 
berg atoms by slow ion-Rydberg collisions |^,|4^, ^9| p3| . 
The same classical approach also explains the anomalous 
autoionization lifetimes of Rydberg electrons in circu- 
larly polarized microwave fields Q and the dynamics 
and stability of circular Rydberg states (i.e., states with 
maximum angular momentum) in weak, slowly rotating 
electric fields |,[50j. Clearly, the success of the classical 
treatment must stem from some special equivalence be- 
tween classical and quantum predictions in the hydrogen 
atom in weak external fields: in fact, such equivalence 
has been already simply observed (but not explained) 
only in the limited case of states with initial zero angular 
momentum [p5, 



In this work we show that the accuracy of the classi- 
cal results does indeed rest on a particularly direct con- 
nection between classical and quantum predictions, and 
demonstrate explicitly that in the perturbative limit the 
quantum expectation values of the angular momentum 
and the Runge-Lenz vector obey exactly the same equa- 
tions as the time averaged classical variables. We also 
investigate the conditions on the fields for which the per- 
turbative approach holds, and we find that classical me- 
chanics seems to require more relaxed constraints on the 
external fields than quantum mechanics, as the classi- 
cal condition for the conservation of the Kepler energy 
-see below- is not equivalent to the quantum condition 
for negligible intermanifold mixing. However, by study- 
ing in detail the contributions to the dynamics which 
stem from n-mixing, we demonstrate that the pertur- 
bative equations of motion for the quantum expectation 
values remain accurate also under the less restrictive clas- 
sical conditions, as long as the dynamics is time aver- 
aged over a Kepler period. In fact, the time averaged 
equations describe well the dynamics of quantum expec- 
tation values also when the Rydberg electron is initially 
excited in a superposition of hydrogenic n-manifolds, i.e., 
its initial state is not a stationary eigenstate of the unper- 
turbed Hamiltonian, but a time dependent wave packet. 
Interestingly, time averaging is precisely the same proce- 
dure which leads to the classical perturbative equations, 
and therefore our result is an explicit example of how 
the scrambling of the principal quantum number (quan- 
tum decoherence) brings about a more direct quantum- 
classical correspondence. Finally, we show that for spe- 
cial superpositions of the eigenstates of the bare Hamil- 
tonian (that is, elliptic states pE, 67,p8|) and in the limit 



tial conditions. Most importantly, elliptic states are not 
merely theoretical constructs: they have been prepared 
in the laboratory and some of their properties have al- 
ready been studied experimentally j|, |69pf7l| . 

Our findings are not merely an application of Ehren- 
fest's theorem Jt^]. Although Ehrenfest's theorem relates 
the time evolution of the quantum expectation values to 
the classical equations of motion, it does not state that 
quantum expectation values and classical variables obey 
exactly the same equations. Such an identity holds only 
for the harmonic oscillator and, albeit trivially, also for 
the free particle and any linear potential. Because the 
harmonic oscillator is often used as a textbook example, 
it often leads to the incorrect impression that such exact 
equivalence is of more general nature. In truth this corre- 
spondence is a very special property of potentials which 
are at most quadratic, because in general the expectation 
value of the "force operator" /(f) = — VV(r), which is 
a function of quantum observables, is not equal to the 
same function evaluated at the expectation values of the 
observables; that is (in general): 



(i) 



In the very special case of the harmonic oscillator, for 
example, the restoring force of the spring is simply pro- 
portional to x and the two sides of Eq. ([!]) are identical, 
hence the exact correspondence between quantum and 
classical evolution. Also, the perturbative treatment of 
the Kepler problem borrows heavily from the methods of 
celestial mechanics |73j and studies the dynamics of time 
averages of the classical variables. Such an approach is 
not the same as an expansion of the Hamiltonian around 
an equilibrium point and up to quadratic terms in the po- 
tential, which would make the system trivially equivalent 
to a harmonic oscillator. Therefore, our work amounts 
to an extension of Ehrenfest theorem, in a much stronger 
form, for the important case of the hydrogen atom in 
weak, slowly varying external fields. 

This paper is organized as follows: in section II we de- 
rive explicitly the equations of motion for the quantum 
expectation values over states confined within a hydro- 
genic n-manifold, and show that they coincide with the 
classical perturbative equations to first order in the fields. 
In section III we investigate the contributions to the dy- 
namics due to the intermanifold mixing: we show that 
the same perturbative equations of motion remain accu- 
rate even if the state is not initially confined within a 
specific n-manifold, as long as one considers the time av- 
erage (over a Kepler period) of the dynamics. In section 
IV we study the initial conditions for the quantum ex- 
pectation values over different quantum states and also 
discuss a few physical implications of our results. Finally 
in section V we draw some general conclusions. 



of large principal quantum numbers, the quantum expec- 
tation values also have the appropriate quasiclassical ini- 
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II. CLASSICAL AND QUANTUM EQUATIONS 

In atomic units (which we use throughout this paper) 
the Hamiltonian for a hydrogen atom in crossed electric 
and magnetic field is: 



H = 



2 



uj l L, 

r 



+ Fx + ^(x 2 + y 2 ), (2) 



where the electric field is parallel to the x-axis and its 
strength is F; the magnetic field is antiparallel to the 
z-axis and ujl is the Larmor frequency of the magnetic 
field, which in atomic units is equal to half the strength of 
the field. For weak fields the diamagnetic term, which is 
proportional to the square of the field, can be neglected, 
and the Hamiltonian becomes identical to the one for 
a hydrogen atom in a weak electric field rotating with 
frequency uil , in the noninertial frame which rotates with 
the field 

The classical perturbative treatment of the hydrogen 
atom in weak, external fields is based on the methods of 
celestial mechanics, and one is interested in the secular 
variation of the elements of the Kepler ellipse followed by 
the classical electron |7j| . While a Kepler ellipse can be 
described by many equivalent sets of elements, the equa- 
tions of motion are particularly simple if one chooses the 
angular momentum and the Runge-Lenz vector. There- 
fore, the dynamical variables of the classical problem are 
the time averages over a Kepler period, and along a Ke- 
pler ellipse (which is the classical solution to zeroth order 
in the external fields), of the angular momentum L of the 
electron, and its scaled Runge-Lenz vector a, which for 
bound states is defined as 1741: 
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^f^2E 12 



(p x L — L x p) 



(3) 



where E = —l/2n 2 is the Kepler energy of the electron. 
The antisymmctrization of the cross product is not nec- 
essary in classical mechanics, but is essential in quantum 
mechanics to obtain a hermitian operator |7^,^6|. To 
first order in the external fields the classical, time aver- 
aged angular momentum and scaled Runge-Lenz vector 
(which for the sake of a simpler notation we will still 
indicate respectively as L and a) satisfy the following 
equations of motion Q fT7|-|8rj : 



dL 

~dt 
da 



= —ll>s xa + ui,xL 



-u>s x L + x a 



(4) 



where ujs is the Stark frequency of the electric field, and 
is defined as: 



-riF 



(5) 



and is the Larmor frequency vector: it is directed 
along the external magnetic field and its magnitude is 
equal to the Larmor frequency of the field itself. 



Eqs. (|i|) were derived originally by Born j77j and in 
classical mechanics they remain accurate as long as the 
two frequencies (Stark and Larmor) are much smaller 
than the Kepler frequency u>k of the electron ^l] : 



UJ S , LU L < L0 K = — 



(0) 



In classical mechanics Eq. (^) means that the elements 
of the Kepler orbit do not vary significantly over a Kepler 
period, so that the classical electron still moves, to a good 
approximation, along a Kepler ellipse, and the Kepler 
energy of the classical electron is conserved. In fact, the 
classical angular momentum and the Runge-Lenz vector 
obey two constraint equations [[74|: 



(7) 



It is easy to see that both constraints are invariant under 
the time evolution dictated by Eqs. (||), and also that 
the second of Eqs. (0) implies the conservation of the 
Kepler energy of the electron. 

Instead, the quantum mechanical interpretation is that 
the external fields remove the degeneracy of the n 2 un- 
perturbed states of the hydrogenic n-manifold, and the 
conditions of Eq. (^) mean that the energy difference 
between two adjacent perturbed states is much smaller 
than the separation between adjacent, unperturbed Ry- 
dberg energy levels. However, this is not the usual con- 
dition under which in quantum mechanics inter- n mixing 
is negligible. For example, in the case of just an external 
d.c. field (the extension to include also a magnetic field is 
straightforward (8^]) the energy separation between the 
lowest and the highest Stark states for a fixed principal 
quantum number n is (to first order in the field) |82j : 



AE = 3n(n-l)F 



(8) 



Therefore the approximate condition for level crossing of 
the highest Stark state from a given n-manifold with the 
lowest Stark level from the next n-manifold is given by 
the Inglis- Teller limit §: 



Zn 2 F 



1 



(9) 



Clearly, in the semiclassical limit the quantum condition 
on the external fields for negligible n-mixing of Eq. (||) is 
much stronger than the classical condition of Eq. (|g) . We 
show below, however, that the perturbative treatment of 
the dynamics of the quantum expectation values remains 
accurate also in presence of some degree of n-mixing in- 
duced by the external fields, as long as the dynamics is 
time averaged over a Kepler period. 

In this section we confine our study to the dynamics of 
the quantum expectation values of the angular momen- 
tum and the scaled Runge-Lenz vector operators (that is, 
(i/}\h\ip) and throughout this paper we use bold- 

face letters for vectors, and a caret indicates a quantum 
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operator, not a unit vector which we denote instead as 
e,) over superpositions of the hydrogen atom eigenfunc- 
tions with a well defined principal quantum number; i.e., 
over states \ip n ) which are localized within a hydrogenic 
n-manifold. More precisely, we show here that to first or- 
der in the external fields the quantum expectation values 
obey exactly the same equations of motion as Eqs. (jij). 

The cornerstone of the study of hydrogenic systems in 
weak, external fields is the so-called Pauli's repacement 
|75| , |7q| , according to which the matrix elements between 
states within the same n-manifold of the position opera- 
tor are directly proportional to the corresponding matrix 
elements of the Runge-Lenz vector operator: 



(nl m \ti\ntm) = — — n(p,t' m!\a l \nlm) 



(10) 



By replacing the position operator, which appears in 
the perturbation Hamiltonian for an external field, with 
— Znk/2 the demonstration of quantum-classical equiva- 
lence within an n-manifold is straightforward |j3f|, |5f|. [B3f , 
but nothing can be said about the intermanifold dynam- 
ics. Instead, in our argument we do not apply Pauli's 
replacement directly in the Hamiltonian; our approach is 
more complicated, but it makes possible the extension of 
our analysis (in the next section) also to the dynamics 
of Rydberg wave packets. Moreover, we will be able to 
show that the quantum analogs of Eqs. ^ remain accu- 
rate under the more relaxed, classical conditions on the 
external fields given in Eq. ([|). 

To prove the special quantum-classical equivalence, we 
will make use of the following identity, which holds within 
a hydrogenic manifold and which we derive explicitly in 
Appendix A: 



field, in a noninertial frame rotating with the field itself 
|p7|,^7|-|50|l). Both the scaled Runge-Lenz vector and the 
angular momentum commute with the hydrogenic Hamil- 
tonian (they are invariants of the pure Kepler problem) 
[ fFij . Moreover, it is easy to see that the magnetic term of 
the classical equations can be recovered by invoking the 
vectorial properties of a and L, because of which their 
commutators with the magnetic term in of the Hamilto- 
nian (i.e., —ijJlLz) obey the well known rule |)72|: 



V u Lj 



(13) 



where V, stands for the i th component of any vector op- 
erator. 

Therefore, we only need to investigate the commuta- 
tors of a and L with the electric perturbation Fx. 
We begin with a y : 



-iF[a y ,x] = -i-T^F { (p z L x - p x L^j ,x 

(tzPx - L x p 7 )j , X I = -^TT^I 



L 7 .x 



Px 



L z ,x 



- [p x ,x] L z - L z [p x ,x] 
= nF \L Z - yp x } . 

However, using the identity of Eq. (|ll]) one has: 

- (lpn\yPx\tl>r, 

from which it follows: 



(14) 



~{ipn\xPy ~ VPx\^n 



{^n\fiPj\lpn) = -(4>n\Pzrj\lpn) 



(11) 



iF(tp n \ [a y ,x] |V>n) = -jnF(ip n \L z \ip n ) 



(15) 



(16) 



where f t and p 3 are components of the position and mo- 
mentum operator respectively. 

Indeed, armed with the result of Eq. |ll| it is easy to 
show that to first order in the external fields the expecta- 
tion values of the quantum observables satisfy the same 
equations as the time averages (over a Kepler period) of 
the classical variables. 

The equations of motion of the quantum expectation 
values are straightforward in the Heisenberg picture p2| : 



(V'nl^rl^n 



W'nl^lV'n) = ~i(lpn\ 



L t ,H 



IV'n 
\lpn, 



(12) 



and we now show that they are identical to Eqs. (|4|). 

The classical equations (|^) contain two terms, an elec- 
tric term which is proportional to the Stark frequency 
uJSi and which couples the angular momentum to the 
scaled Runge-Lenz vector; and a magnetic term which 
is proportional to the Larmor frequency of the field 
(or to the rotation frequency of a slowly rotating electric 



This is the same as the electric term in the equation of 
motion for the classical time averaged a y . 

The derivation of the electric term for the equation of 
of {ip n \o- z \^n) follows along the same lines and it is easy 
to see that it yields the desired result. 

We then consider a x : 

-iF[a x ,x] = (p y L z -p z Ly^ ,x 

\LyP Z ~ L Z Py) ,xj = -l^F \f>y 



L z , x 



Ly 7 X 



Ly, x 



L-,,x 



Py 



) 



(17) 



2 F \PyV + VPy + P* z + z Pz s 



Invoking once again Eq. ( |11| ) one obtains immediately: 

F<^„|[ax,x]|^„) = 0, (18) 

which is the same as the right hand side of the corre- 
sponding classical equation of motion. 
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Finally, we turn to the equations for the angular mo- 
mentum. The classical equations can be written as: 



dt 



-e ltJ tws.Ofc = ——nFet i kCLk 



(19) 



where we have specialized the right hand side to the case 
of an external field along the x-axis. Using once again 
the vector properties of f , the contribution of the electric 
term to the quantum equations is: 



iF L 



Fe, llk fk 



(20) 



This is not yet in the desired form. However, within 
a given n-manifold one can apply Pauli's replacement 
|7q,[7q] according to which: 



(ni'm'lfilnim) — — nijil' rri\a l \nlrri) 



(21) 



Pauli's replacement is mathematically exact, and yet 
physically it is just an approxximation, because the dy- 
namics of the electron is only approximately confined 
within a given n-manifold. Clearly, the accuracy of the 
approximation rests on certain conditions on the exter- 
nal fields, depending on which the dynamics may -or may 
not- be very well localized within a hydrogenic manifold, 
and we discuss such conditions in detail below. However, 
in the present section we are interested only the intra- 
manifold dynamics, and therefore one has: 

Fe lA}k (fpn\fk\ipn) = -7;nFe lA:k (ip n \ak\ipn) , 

2 (22) 

which is the same electric term as in the right hand side 
of the corresponding classical equations of motion, and 
our proof is complete. 

Our derivation of the equations of motion of the quan- 
tum expectation values is accurate only to first order in 
the fields because it relies heavily on the identity of Eq. 
(|TT|), which holds for the unperturbed \tp n ) states; these 
are eigenstates of the hydrogen atom Hamiltonian and 
therefore are quantum solutions to zeroth order in the 
fields. In fact, this is exactly in the same spirit as the 
classical approach, where the right hand sides of Hamil- 
ton equations are time averaged over a Kepler period 
and, most importantly, along Kepler ellipses, which are 
the classical solutions to zeroth order in the fields. 

The same consideration can also be cast in the lan- 
guage of operators, by observing that we have proven 
Eq. ([ll]) only for the time independent operators of the 
Schrodinger picture, whereas in the right hand side of the 
equations of motion one must more correctly use the time 
dependent operators of the Heisenberg picture. However, 
in a first order approximation one may assume that the 
time evolution of the operators in the Heisenberg picture 
is dictated solely by the hydrogenic propagator, which 
commutes with both L and a, so that one can legitimately 
use the properties of those operators in the Schrodinger 
picture. 



Since Eq. (|11|) holds for all possible pairs of indexes 
{j, k}, the same derivation can be easily extended to the 
case of slowly varying (both in magnitude and direction) 
electric and magnetic fields, in which case the perturbing 
Hamiltonian H± becomes: 



H 1 = '%2{F % (t)t t -u Lt (t)L % } 



(23) 



However, the most important feature of our proof is 
that we have not applied Pauli's replacement directly 
in the Hamiltonian pi], [51, 83 1. An early application 
of Pauli's replacement yields a straightforward proof of 
quantum-classical equivalence within an n-manifold, but 
it erases all information about the precise conditions on 
the fields under which the perturbative classical equa- 
tions of motion constitute an accurate description of the 
dynamics of the quantum expectation values. Moreover, 
it also makes impossible to study the corrections to the 
dynamics due to intermanifold mixing and therefore one 
could not extend Eqs. (|j) to the case of Rydberg wave 
packets. Instead, in the next section we address in detail 
precisely these important issues. 



III. INTERMANIFOLD DYNAMICS AND 
QUANTUM-CLASSICAL CORRESPONDENCE 
IN RYDBERG WAVE PACKETS 

In this section we study intermanifold mixing and the 
conditions on the external fields under which the clas- 
sical perturbative equations of motion offer an accurate 
treatment of the dynamics of the quantum expectation 
values. We will show that the same conditions as in clas- 
sical mechanics hold in quantum mechanics too, provided 
that the dynamics is time averaged over a Kepler period. 
Most importantly, we will demonstrate that upon time 
averaging Eqs. (||) apply also to the case of Rydberg 
wave packets. 

For the sake of simplicity we restrict our analysis to the 
pure Stark case, that is when there is no external mag- 
netic field; the extension to the more general case includ- 
ing a weak magnetic field is straightforward. Therefore, 
in this section we assume a simplified Hamiltonian: 



H = H Q + Fx 



(24) 



where Hq is the hydrogen atom Hamiltonian and F is 
again the strength of the external electric field. 

Inter-manifold mixing is due to two main causes, de- 
pending on how the Rydberg state is prepared. If the 
Rydberg electron is initially confined within a hydrogenic 
n-manifold, then n-mixing is induced by the applied ex- 
ternal field, and in that case the intermanifold contri- 
butions to the equations of motion are of second order 
in the applied field. Note that this is a very realistic 
picture for slow ion- Rydberg collisions ^,^8|,|5^ In 
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fact, ion-Rydberg collisions are actually gentle encoun- 
ters at very large ion-Rydberg separation, which are ef- 
fective because of the long-range nature of the coulomb 
interaction, and are very accurately modeled by a time 
dependent, weak external field acting on the Rydberg 
electron |47| ) |4q] . Typically, the Rydberg state is pre- 
pared in absence of external dc fields, and the weak field 
of the colliding ion is turned on adiabatically as the ion 
slowly approaches the Rydberg atom. In such situation, 
the analysis of intermanifold mixing is equivalent to the 
study of the second order corrections to the equations of 
motion, and it allows one to determine for what precise 
conditions on the external electric field such corrections 
are negligible. 

However, n-mixing can also be present at the outset, 
either if the Rydberg state is prepared in presence of the 
applied dc field and the field is strong enough to mix ad- 
jacent Rydberg levels, or alternatively if a short, large 
bandwidth laser pulse is employed in the preparation of 
the Rydberg state. In both cases the ground state is 
coupled to a distribution of hydrogenic manifolds, and 
the Rydberg electron is not excited to a high energy, sta- 
tionary eigenstate of the hydrogen atom, but rather to 
some time dependent wave packet. To a first approxima- 
tion the Rydberg wave packet oscillates with the Kepler 
frequency of the eigenstate around which the distribu- 
tion of principal quantum numbers is centered; therefore 
a Rydberg wave packet contributes rapidly oscillating, 
intermanifold terms to the equations of motion for the 
quantum expectation values. Such intermanifold contri- 
butions may be of first order in the external field. In 
this section we extend the validity of the classical, per- 
turbative Eqs. (Q) precisely to the case of Rydberg wave 
packets, by time averaging the equations of motion over 
a Kepler period and by showing that the secular, inter- 
manifold contributions to the dynamics remain negligible 
under the classical conditions of Eq. (JsJ) for the external 
field. 

More precisely, the wave function of a Rydberg wave 
packet is: 

Mi)) =£C- n |Vn(*)} . (25) 

n 

and the Heisenberg equations of motion for the expecta- 
tion values of either L or a over the state of Eq. ( pq ) 
include "off-diagonal" matrix elements (for the sake of 
brevity we call "off-diagonal" the matrix elements of any 
operator between states from two different hydrogenic 
manifolds; whereas we will call "diagonal" the matrix ele- 
ments between two states within the same n-manifold, re- 
gardless of their angular momentum quantum numbers). 

For example, if we indicate generically by O c i 4- O q the 
combinations of operators in the right hand side of the 
Heisenberg equations of motion, in the case of L r one has: 

j t (tJj\L t (t)\i>) =F[J2 \°n\ 2 (i>n\Ocl(t) + O q {t)\^n) 



+ Y, Cn>c n {ii> n >\d cl {t) + 6 q {t)\iij n )} . 

n ,n 

where O c i indicates a combination of operators which 
correspond to the classical variables in the right hand side 
of the classical equations of motion (i.e., a combination of 
angular momentum and Runge Lenz vector). Instead, O q 
indicates the purely quantum corrections, which vanish 
when the motion is exactly confined within a hydrogenic 
manifold. Obviously, an identical expression holds also 
for the equation of motion of the expectation value of a x 

The expectation values of the operators of O c i appear 
in the equations of motion in the same way as the clas- 
sical time averaged variables, and therefore they evolve 
in time exactly like their classical counterparts. This is 
true even if the state is not confined within a hydrogenic 
manifold. Yet, for states which are spread over more than 
one n-manifold the quantum contributions from the ex- 
pectation value of O q do not vanish exactly. In fact, we 
proved in the previous section and in Appendix A that 
only the diagonal matrix elements of O q vanish, but the 
same does not hold for the off-diagonal matrix elements. 

However, we show below that under the classical con- 
ditions for the external field, and upon time averaging, all 
the off-diagonal matrix elements of the double sum of Eq. 
( p6| ) offer a negligible contribution to the dynamics. In 
our demonstration we do not distinguish between the two 
operators, and treat O c i + Oq as a single term which, for 
the sake of brevity, we simply denote as O. That is, our 
argument demonstrates that also the off-diagonal matrix 
elements of O c i yield only negligible contributions to the 
equations of motion. Therefore, the quantum dynamics 
of the expectation values of the operators of O c i is deter- 
mined only by the intramanifold terms. Note that this 
last observation is not essential to the issue of quantum- 
classical correspondence. 

Before our demonstration, however, we must discuss 
briefly the operators of Eq. (|2^) and the magnitude of 
their matrix elements between two eigenstates of the hy- 
drogen atom. From the previous section, and also from 
Appendix A, it is easy to see that the sum of the two 
combination of operators O c i + O q = O may be equal to 
one or a combination of the following operators: 

O ~ i np % f 3 i, j = 1,2,3 , (27) 
[nL l 

First, the nL l operator yields matrix elements the magni- 
tudes of which are at most < n 2 . Moreover, L l commutes 
with the hydrogen atom Hamiltonian and therefore all its 
off-diagonal matrix elements vanish, and it is very easy 
to prove that its intermanifold contributions to the dy- 
namics are negligible (see below). 

Next, in the case of f, the magnitude of the matrix 
elements and, most importantly, their scaling with n are 



G 



determined solely by the radial matrix elements: 

Rl f = (nlm\r\n'l'm') . (28) 
If n' = n one has 182]: 



R 



ni- 



_ r>n,t.- 
1 — ^nJ 



1 = --n 2 
2 



1 



f 2 



(29) 



whereas if n' ^ n the radial matrix elements of the 
position operator are given by a complicated formula 
which involves hypergeometric functions [p2|. However 



for n', n 3> 1, 



is accurately approximated by a 



well known semiclassical result 



Rl 



1 + Ao' — MfA (A n „<e) 
n,- ' 



(30) 



< 



2A, 



where A n! „/ = n — n' and A^/< =£ — £'. In Eq. (30) we 
also used: n c = 2nn' /(n + n 1 ), l c — ma,x(£,C) and finally 
e 2 = 1 - t 2 c /n 2 . Clearly the last inequality of Eq. ( |30| ) is 
accurate only to the leading order in n, and to the same 
order it is equally correct if one uses either n or n' or n c . 



and insert- 



Finally, for np^j by using p t — — 

ing a resolution of unity between the two operators one 
has: 

\n (n'e'm'\p l f 3 \n"£"m")\ 



in y (n'£'m'\f t \fi) (n\?,\ri'£"m") 



< \ (n'£'m'\r J \n"e"m")\ 



< 



(31) 



where A n /^ — n' — n u , and where we used the final result 



of Eq. (|3C|) . In Eq. (|3lJ) for the sake of a simpler notation 
we adopted the following convention: \fi) = {n^^m^} 
(which we will often use in this section). 

The result of Eq. (|3l| ) rests on the observation that 
in the semiclassical limit the radial matrix elements of r l 
become rapidly very small for large A's, as one can easily 
see from Eq. (|30"|); and one may safely assume that for 
nonnegligible matrix elements the difference between the 
two principal quantum numbers is always much smaller 
than any of the principal quantum numbers themselves 
(e.g. A n / !jU <^ n',n^). Therefore one may legitimately 
neglect higher order corrections in A„/ jM /n'. 

In fact, we assume precisely this important condition 
throughout our argument: that is, we assume that the 
variance of the distribution of the Rydberg wave pack- 
ets over the hydrogenic principal quantum numbers is 
always much smaller than the principal quantum num- 
ber at the center of the distribution, i.e., the approximate 
average principal quantum number of the Rydberg wave 



packet. This is a very realistic approximation for most 
laser pulses employed in the excitation of Rydberg elec- 
trons, and it breaks down only for ultrashort, ultralarge 
bandwidth pulses; or when the Rydberg state is excited 
in presence of ultrastrong external fields. 

In what follows we conduct our analysis in the most 
general form. Although at some point we specialize our 
argument to the case O = r l which yields the largest off- 
diagonal contribution to the equations of motion, it will 
be easy to see that the treatment of the case O = np % r 3 
is completely analogous. 

Any off-diagonal matrix element of Eq. (|2^) can be 
written as: 



(iP n ,\e iAt de- iAt \ip n ) = ^ £ ,m')C n (l,m) 
x (n'e'm'\e i6t 6e- i6t \nem) , ^ 



where we have expanded the states |Vv) and \ip n ), which 
are initially confined within the n'- and ?i-manifold re- 
spectively, as follows: 



Hn) =^C n {i,m)\nlm) 



(33) 



where the C n (£, m)'s are some general coefficients, possi- 
bly complex. 

As we mentioned before, these off-diagonal matrix ele- 
ments are present in the equations of motion either when 
the high-energy electron is prepared in a wave packet, or 
when the electron is initially confined within a single n- 
manifold, in which case they represent the second order 
(in the external field) corrections to the dynamics. 

The dynamics of the matrix elements of Eq. (^) is 
best studied in the interaction picture J72] | , and therefore 
it is convenient to set: 



e- i6t \n£m) = 53 ax{t) e - lE ^\X) 

A 



(34) 



where we used once again the convention |A) = \n\£\m\) 
and = In^i^rrifi), as the most important features of 
our argument depend on the spectrum of the hydrogen 
atom, and are determined solely by the principal quan- 
tum number of the state. 

The equations of motion for the a's and the /?'s can be 
derived directly from the Schrodinger equation: 



tax 



= Fj2(Mx\M)ax ie - 



i(E Xl -E x )t 



(35) 



and the solution of Eqs. (|3J]) to zeroth order in the field 
is: 



a 



x ,l^m x ,m 



(36) 
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In a first order approximation Eqs. ( |35| ) become: 

ia[ F) = F(X\x\nem)e-^ E "- E ^ t 



ipl F ) = -F(n'£'m'\x\n)e 



(37) 



and their solution is straightforward: 

a f ] = F(X\x\n£m) e 'g^,^' -1 n x ^ n 

— —iF(\\x\n£m)t ri\ — n, {^a ; to a} 7^ {(■, m } 
A F) = 1 {n A ,^,m A } = {n,^,m}( 38 ) 



and also: 



^ = F{n'i'm'\x W -^f^g^ 



^ n' 



fiP = iF{n'£'m'\x\fx)t n„ = n', {^,m M } + {£', m'j 



Therefore, to first order in the external field the time 
dependence of the matrix elements of Eq. (32) is: 



stF) 



(n'£'m'\e lNt Oe- lHt \n£m) 
= Fi^(nYm'|0|A)(A|x|n^ 



-i{E n -E n ,)t 



-i{E x -E„,)t 



X 



E n — E\ 
F 2 Y}n'£'m'\x\^) ((j,\6\nlm) 

p -i(E n -E n ,)t _ e -i(E n -E„)t 

E n i — E„. 



— iFs'^^(n'£'m'\d\n£\m\) (n£\m\\x\n£m) 

x e -i(E n -E n ,) H 

+ FiY^{ n 'Z' m '\6\n'l\m\) (n'£ x mx\x\n£m) 

e -i(E n -E n ,)t _ 1 



(40) 



X 



>iFz^{n'i'm'\x\n'£»mp){n'l»mn\6\nlrn) 

e -i(E n -E n ,)t t 

Fe'y^^n'£'m'\x\n£ ll m fl ) (n£ ^\0\n£m) 

e -i{E n -E n ,)t _ 1 



En 



+ (n'£'m'\6\n£m) 



^-i{E n ~E n ,)t 



where we have attached subscripts to the field strength 
F only for bookkeeping purposes, and so F± = F% = ... = 
F 6 = F. 

Clearly, all the terms of Eq. ([hj|) oscillate with a fre- 
quency comparable (but not identical!) to the Kepler 



frequency ujk of the hydrogenic manifold at the center 
of the distribution of principal quantum numbers. We 
indicate the principal quantum number of this special 
hydrogenic eigenmanifold as n. For very weak external 
fields, the Kepler frequency is much larger than the Stark 
frequency of the motion. This is true for the classical 
conditions on the fields of Eq. (^) and also, albeit in 
a much stronger form, for the usual quantum condition 
of Eq. (^), i.e., the Inglis- Teller limit for negligible in- 
termanifold mixing. This means that the wave packet 
oscillates several times before the classical perturbative 
equations of motion yield any significant change in the ex- 
pectation values of angular momentum and Runge-Lenz 
vector. Therefore, following exactly the approach of clas- 
sical perturbation theory for the derivation of Eqs. (Q), 
we time average the quantum dynamics over the Kepler 
period Tk = litv? of the n-manifold approximately at 
the center of the energy distribution of the wave packet. 
Such time averaging does not affect the diagonal terms of 
Eq (|2^); however, it allows us to evaluate with accuracy 
the secular, off-diagonal contributions to the dynamics 
over a time r = 7T5, that is, over some multiple of the 
Stark period T5. 

For a Rydberg electron confined within a n-manifold 
the Stark period Ts is defined as: 

Ts = 3^ ■ < 41 > 
However, its definition can be easily generalized to the 
case of a wave packet with average principal quantum 
number n, by setting: 

Ts ~ ir ■ (42) 



First, we consider the last term of Eq. (40) which is 
of zeroth order in the field; by inserting it in the pertur- 
bative equations of motion, i.e., in Eq. (p6|), one obtains 
an off-diagonal term which is of first order in the field. 
By first time averaging over a Kepler period, and then 
integrating over a time r, and finally using Eq. (^2] one 
obtains: 



F{n'\6\n) 



e ^~ E ^\ dr' 

K 



« l^=(n'\6\n)^ (a^ ~ ^A n ,^ , (43) 

where (■■■) K indicates the time averaging over the Kepler 
period, and where we used the following result (derived 
in Appendix B): 



-i(E n -E n ,)t 



1 

K ~ Tk Jo 



,Ti(E,-E 3 )t dt 



-A 



O 



A- 



From Eq. (pp|) one has: 



\{n'\d\n)\ = \{n%\n)\ 



< 



(44) 



(45) 
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which holds only to the leading order in n, and where a 
similar result can be obtained also in the case O = np % rj 
The final result is: 



< 



i(E„-E n ,)t\ 

2A^_ n A n / n 



F (n'lOln)^-^"-^) dr' 



2|A„, J 



(46) 



< n , 



where the last inequality follows under the very impor- 
tant assumption which we introduced before, i.e., that 
the variance of the distribution of the Rydberg wave 
packet over the principal quantum numbers is much 
smaller than the average principal quantum number of 
the wave packet. We conclude from Eq. (|4^) that first or- 
der, off-diagonal contributions are very small when com- 
pared to the dominating principal quantum number of 
the wave packet; indeed, this is a sufficient condition to 
neglect them completely, because the quantum expec- 
tation values of angular momentum and scaled Runge- 
Lenz vector range precisely from —n to n. Most impor- 
tantly, however, one must also require that 7 ~ 1, i.e., 
the off-diagonal contribution remains small only up to 
times comparable to the Stark period; for longer times 
the secular effects build up and off-diagonal terms be- 
come relevant. 

Clearly, the analysis above does not yet yield any in- 
formation about the precise condition that the external 
field F must satisfy, so that all off-diagonal terms remain 
negligible. To learn more about it, one must analyze the 
off-diagonal contributions which arc of second order in 
the field. We begin considering the terms of Eq. ( p0| ) 
which are proportional to F±. By inserting any of them 
in the equations of motion for the expectation values of 
angular momentum and Runge-Lenz vector, one obtains 
contributions which are of second order in the external 
field; that is, to the leading order in n one obtains: 

( e -i(E n -E n ,)t\ 

F 2 (n'£'m'\6\X)(X\x\n£m) dr' 

! 1 E n — E\ 



< 



A* 



,<n> 



2A r 



1 n.A 



Fn b 



(47) 



where Jn) = (n + n')/2, and we have used the results of 
Eqs. (p2|), ( p4| ) and (|4|); we have also used the following 
result (see Appendix B): 



1 



E, 



« U 3 ( \ 1 . 



-0 



A^ 

ft 2 



(48) 



From Eq. (|47]) is finally possible to extract a necessary 
and sufficient condition on the field strength. By requir- 
ing that the result of Eq. (^) is much smaller than n 
one obtains: 



which is essentially the same as the classical condition of 
Eq. (^), as we had claimed before. 

Clearly, the same analysis applies to the terms of Eq. 
( [f0| ) which are proportional to F2, and also to the oscil- 
lating part of the F4 and Fq terms. 

Therefore, we next turn our attention to the remaining 
terms from Eq. (40), that is, the contributions which are 
proportional to F3 and F§ and also the non-oscillating 
parts of the F4 and F$ terms as well. First, one needs 
the time averages of the time dependent factors in the F3 
and -F5 sums, which are given by the following equation 
(see Appendix B): 

it e -i&-E,)t\ = — ( TK t e-^-^dt 



>K Tkj 

n 3 ( 6 / f 
-— { 1 - - ( An, 3 - ~A„ ) (f + irXjj 



-O 



A 3 



/1 



(50) 



By inserting the time averaged F3 and .F5 terms, along 
with the non oscillating parts of the F4 and Fq terms in 
the equations of motion and integrating over time one 
obtains (to the leading order in n), an additional inter- 
manifold correction, which we denote as G(n'£'m';n£m): 

-3 

G(n'£'m'; n£m) = 7 T S F 



F 3 (n'£'m'\d\n£xmx)(n£ x mx\x\n£m) 
F A ^2 (n'£'m'\d\n'£ x m x )(n'£ x mx\x\n£m) 
F 5 Y (n'£'m'\x\n , £ f ,Tn fl )(n'£ tl m IM \d\rdm) 



+F 6 (n't'm'lxlnt^m^int^m^Olrdm) 

ft is easy to see from the previous analysis of the F\ 
terms that G(n'£'m';n£m) is a negligible intermanifold 
contribution if and only if the expression within the curly 
brackets of Eq. (|5(]) scales as ~ Fn 3 . However, each of 
the terms in the four sums consists of the product of 
two matrix elements, and our previous analysis of the 
magnitude of the matrix elements of r l and also of any of 
the possible choices for O, indicates that all such terms 
may scale as ~ Fn A . Therefore, the desired scaling as the 
third power of n must originate from cross cancelations 
between the four sums within the curly brackets. 

This is obviously correct when, for example, O = 
nL z and all the off-diagonal matrix elements of 



Fn 1 < f 



(49) 







G(n'£'m';n£m) vanish, and one has: 

\G(n'£'m';n£m)\ 
2n 2 

= 7tt-t \\F n ( m ~ rn'){n'£'m'\x\n£m}\ 
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Fn 5 



(52) 



where we used Eq. (|42), and also the usual selection rules 
to conclude that m — m! = ±1. A moment's thought 
shows that an essentially similar analysis holds also for 
the other components of the angular momentum. 

The situation, instead, becomes much more compli- 
cated when O = r l or O = np^fj In those cases the 
pairing of terms which leads to the desired cross cancela- 
tions depends on the differences of the angular momen- 
tum quantum numbers; that is, it depends on A^ t e and 
possibly also on A m >. m . For example, if £ = £' — 2, by 
pairing an F3 term with the corresponding F4 term one 
obtains: 



where 77 is a coefficient ~ 1 which contains the angular 
part of the matrix elements |g2[ . The difference of the two 
diagonal radial matrix elements scales as ~ n ~ n, and 
therefore the whole expression scales as ~ Fn 3 , which is 
the desired result. Note that in Eq. ([53|) one can approx- 
imately factor out an off-diagonal radial matrix clement 
because both off-diagonal elements represent the same 
kind of transition, i.e., the principal quantum number 
and the angular momentum quantum number decrease 
in both cases (assuming that n' > n), that is the transi- 
tions are: 



E ■ (53) 



£'^£ x £'>£ x 
h -> t £x>£. 



(54) 



Indeed, it is well known that the matrix element for an 
atomic transition which increases both the energy and 
the angular momentum of the electron is significantly 
larger than the one for a transition which brings about 
the same change in energy, but leads to a smaller final an- 
gular momentum f&q] . Therefore, if A^^ = 0a different, 
more complicated pairing of the terms must be employed, 
which in this case may depend also on the angular part 
of the matrix elements. 

However, the scaling of G(n'£'m';n£m) with the prin- 
cipal quantum number can most effectively and also more 
convincingly be studied by evaluating numerically the 
whole expression within the curly brackets of Eq. (51), 
divided by F. More precisely, for each pair {n',n} we 
computed the maximum magnitude of the expression 
within curly brackets (divided by F) over all possible 



choices of angular quantum numbers, and we we denote 
it by g(n', n), that is: 



g(n',n) 



max 

{l'm';tm} 



,G(n'l'm';n£m) 



jTsF 2 n 3 



(55) 
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n 

FIG. 1. Scaling of the intermanifold contribution with the 
principal quantum number. In Fig (a) we plot the natural 
logarithm of the maximum magnitude of the intermanifold 
terms vs. the natural logarithm of the principal quantum 
number, and the approximate straight line indicates a simple 
power-law scaling. In Fig (b) we plot the slope of the line in 
Fig. (a), i.e., the exponent of the power-law, which is clearly 
converging to f = 3, thereby proving that in the semiclassical 
limit all intermanifold contributions become negligible. 

In Fig. |l|- (a) we plot the logarithm of g(n',n' — 1) vs. 
the logarithm of n', setting O = x: this is the case in 
which G(n'£'m'; n£m) is the largest, as many numerical 
calculations confirm. The line in Fig. [j]-(a) is almost 
exactly a straight, which confirms that g{n' , n'—l) scales 
with n' i=a n according to a power law, i.e. g ~ rfi. 
However, the exponent £ is not exactly constant, and in 
Fig. 0-(b) we plot the slope of the straight line of Fig. 
|l|-(a) vs. n' . Clearly, for increasing n"s the exponent £ is 
converging to 3, and that is precisely the result which we 
need to prove that G(n'£'m'; n£m) is truly a negligible 
intermanifold contribution to the equations of motion. 

The proportionality coefficient of the power law can 
also be easily evaluated from the numerical data and one 
obtains: 



g(n', n' - 1) w 3.15n' 3 = 3.15n 3 j 1 + O (^j j 



(56) 



Finally, we have repeated the same calculations for sev- 
eral values of A n / „, and also for O = y, z, and O = np^j, 
{1,]} = 1,2,3. In all cases our findings were essentially 
identical to the ones of Fig. |l| 
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The physical interpretation of our result is particularly 
interesting. Rydberg wave packets which are relatively 
well localized in energy move along the trajectories of the 
classical electron or, for radial wave packets, of an ensem- 
ble of classical electrons [p0|-|43|]. However, unless some 
suitable external fields are applied to the system, |32] - 
p3| , they do so only for a few Kepler periods. Such wave 
packets do not remain localized in the angular variable 
and therefore spread along the classical ellipse. Even- 
tually they display interference fringes as the front of 
the packet catches up with its tail, and finally they also 
show quantum revivals and superrevivals pC 22, 25|,|27| 
p9fl . However, from the point of view of the time aver- 
aged equations of motion the electron is always spread 
(i.e., averaged) along the classical trajectory, very much 
like in classical mechanics, where after time averaging the 
elements of the Kepler ellipse become the dynamical vari- 
ables of the system, replacing the phase space coordinates 
of the electron. Therefore the time averaged, quantum 
equations of motion are insensitive to the spreading of 
the wave packet, and to its revivals, and that is why the 
same results as for stationary states can be legitimately 
extended to wave packets too. 

Note that the results of our analysis can be gener- 
alized beyond the case of Rydberg electrons excited to 
wave packets, and hold also when the Rydberg electron 
is initially confined within a hydrogenic manifold, and 
n-mixing is brought about only by external field. In 
that case, in fact, one needs only to replace (n'£'m'\ with 
(n£'m'\. Clearly, for a Rydberg electron initially con- 
fined within a n manifold, all off-diagonal terms in the 
equations of motion for the quantum expectation values 
are of second order in the field, however they still derive 
from the first order terms of Eq. (|io|), and that is why 
the classical constraint on the external field of Eq. p9| ) 
holds in that case too. 

Indeed, in the very important case of slow ion-Rydberg 
collisions the Rydberg electron is initially excited to a 
specific n-manifold, and in the next section we show 
that some special linear combinations of hydrogen atom 
eigenf unctions which are confined within a hydrogenic 
manifold (elliptic states @^,|||) yield the appropri- 
ate quasiclassical initial conditions for the quantum ex- 
pectation values of angular momentum and Runge-Lenz 
vector, which then closely track the time averages of the 
classical variables. 



IV. CLASSICAL, QUASICLASSICAL AND 
QUANTUM INITIAL CONDITIONS 

In this section we discuss the initial conditions for the 
equations of motion for the quantum expectation values, 
and show that in the case of elliptic states [[29[ |67|, |6§|l 
the quantum expectation values track exactly the time 
averages of the classical variables. More precisely, for a 
classical Kepler ellipse, initially in the a;?/-plane, and with 
the semimajor axis pointing along the z-axis, one has: 



L , 



L„ = 



a x 
a. 



ne 
, 



(57) 



where e is the eccentricity of the orbit, and, as usual, the 
energy of the classical electron is: 



E = 



1 

2^2 



(58) 



On the other hand, for an elliptic state \na) (an elliptic 
state is given by a complicated superposition of spheri- 
cal eigenstates of the hydrogen atom, all with the same 
principal quantum number (2^, ^,^38|), which is also lo- 
calized in the xy-plane, and oriented like the classical 
ellipse above, one has f|||67|,|8): 



(na\L z \na) = (n — 1) cos a 
(na\a x \na) = (n — l)sina 

(na\L x \na) = (na\L y \na) = 
(na\a y \na) = (na\a z \na) = , 



(59) 



so that the correspondence between e and sin a is estab- 
lished (obviously this a has no relation with the coeffi- 
cients of the previous section). Clearly, in the limit of 
large ra's the quantum expectation values and the classi- 
cal predictions converge, and for elliptic states the quan- 
tum expectation values not only obey the same pertur- 
bative equations as the time averages of the classical 
variables, but also have almost the same initial condi- 
tions. Therefore they closely follow the same trajecto- 
ries as the time averaged classical angular momentum 
and Runge-Lenz vector. This result has already been 
observed numerically p , 51 , B3|, and also experimentally 
for some special configurations of the external fields. 
Most importantly, since elliptic states are coherent states 
of the angular momentum p9],|6^,|68|,|85||, i.e., states of 
minimum uncertainty, it turns out that: 



(na\L 2 \na) — ^J^nalL^na) 2 
\ {{na\L 2 \na) + J2 l ( na \Li\na) 2 ^ 



(60) 



(n — 1) cos 2 a + 



and in the semiclassical limit the expectation values of 
ij, i = l,2, 3, are related to the expectation value of L 2 
approximately like in classical mechanics. In fact, it has 
been verified numerically in a few special cases J5l|,|8^] 
that not only the expectation values of L z and a t evolve in 
time quasiclassically according to Eqs. (jj), but also that 
during the time evolution the state remains elliptic. This 
is exactly the same situation as in classical mechanics, 
where the electron keeps moving along an ellipse, but the 
properties of the ellipse vary slowly in time; similarly, the 
numerical evidence shows that the elliptic state remains 
localized along a classical Kepler ellipse while it slowly 
evolves in time. 
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Instead, for the more familiar spherical eigenstates 
\n£m)'s of the hydrogen atom the situation is completely 
different, and one has: 



(n£m\L z \n£m) = m 
(n£m\L x \n£m) = (n£m\L y \n£m) = 
(n£m\a x \n£m) = (n£m\a y \n£m) = 
= (n£m\a z \n£m) = , 



(61) 



and the initial conditions differ dramatically from the 
classical ones, which leads to some interesting considera- 
tions. 

As we explained before, the classical constraints of Eq. 
(Q) remain invariant under the evolution of the perturba- 
tive equations. The second classical constraint translates 
into a condition over the quantum expectation values, 
which is also invariant: 



X 



E{<</^> : 



(62) 



where \ip) is any state (elliptic, spherical or also a wave 
packet, in which case the expectation value must also 
be averaged over a Kepler period) which satisfies the re- 
quirements of our derivation. The invariant x is related 
to the Casimir operator of the SO (4) symmetry group 
of the hydrogen atom [[76], ^6j. The value of the classi- 
cal invariant of Eq. (Q) is nr. For an elliptic state the 
quantum invariant of Eq. ( |62| ) is equal to (n — l) 2 ; how- 
ever, for a spherical eigenstate x is equal to m , which 
for the small m's typically excited by optical transitions 
from the initial low-n state to the Rydberg high-n state 
is a much smaller number than the classical result. This 
poses severe limits on the largest possible expectation 
value of any component of the angular momentum over 
a spherical eigenstate; this feature might be exploited in 
experiments to study the properties of Rydberg states. 

The limitation which Eq. (^2|) imposes on the expecta- 
tion values of angular momentum and Runge-Lenz vector 
over spherical eigenstates stems from the fact that such 
states have vanishing electric dipole moment; a small an- 
gular momentum is not balanced by a large Runge-Lenz 
vector, as it happens in classical ellipses and quantum 
elliptic states. More precisely, for spherical states, the 
electric field cannot induce first order dynamical effects 
because the expectation value of the Runge-Lenz vector 
over a spherical eigenstate is zero (i.e., there is no perma- 
nent electric dipole) , and the state must first be distorted 
by the field so that the expectation value of the angu- 
lar momentum (or Runge-Lenz vector) can change. This 
indicates that the dynamics must be at least of second 
order in the external fields. This situation is germane to 
the well known linear Stark effect (8^] , where degenerate 
perturbation theory and parabolic states must be used 
to account for the linear dependence of the eigenvalues 
on the external field. In fact, a spherical eigenstate can 
be seen as a superposition of elliptic states (or in a semi- 
classical interpretation an ensemble of Kepler ellipses), 



which are oriented uniformly in the i?/-plane, so that the 
total Runge-Lenz vector is averaged to zero. 

The (at least) quadratic dependence of the time evolu- 
tion on the field can be explicitly verified by expanding 
the time dependent operators in the Heisenberg picture, 
|f72| and by showing that the expectation values of the 
first order terms in the electric field vanish. Using the 
Hamiltonian of Eq. (Q) (minus the diamagnetic term), 
and writing the hydrogen atom Hamiltonian as Hq one 
has: 



e iHt L,e- iHt =L-it 



Hq — ujlL z , L l 



- itF 



x. L, 



[h - ui l L z , 



k=0 

L. 



Hq - U) L L Z 



(63) 



0(F 2 ) 



The simplest first order term is a direct commutator of 
the field with a component of the angular momentum: 
therefore it is either zero or a component of the position 
operator (depending on the index i) , and its expectation 
value over a spherical eigenstate vanishes. Other, more 
complex first order terms come from the double sum in 
Eq. (|6^) and consist of a first series of k commutators 
of L % with Hq — ujlL Zi and then of a commutator with 
the field, and finally of a second series of (p — 1 — k) 
commutators with Hq — ulL z . It is easy to see that the 
expectation value of any of such terms over a spherical 
eigenstate of the hydrogen atom which is quantized along 
the 2-axis vanishes, simply because such state is an eigen- 
states of Hq — LO]_,L z - However, the same result can be 
proven as follows also when the atom is quantized along 
an arbitrary direction. I) The first series of k commuta- 
tors either vanishes or yields a component of the angular 
momentum (depending on the index i); II) the commu- 
tator of the result of step I with the electric field either 
vanishes or yields a component of the position operator; 
III) finally, the result of steps I and II must be commuted 
(p — 1 — k) times with with Hq — u>lL z ; this sequence of 
commutators can be organized so that one does first the 
commutators with L z and next those with Hq because 
these two operators commute with one another. Obvi- 
ously, the expectation value of the commutator of any 
operator with Hq vanishes, if it is taken over an eigen- 
state of Hq itself. On the other hand the commutators 
with L z either vanish directly or yield a component of the 
position operator, whose expectation value over a spher- 
ical state also vanishes, and our point is proved. 

Similar considerations apply also to L 2 . In fact, it is 
even easier to demonstrate that the time evolution of the 
expectation value of L 2 over a spherical eigenstate of the 
hydrogen atom is of second order in the field. Clearly, 
\n£m) is an eigenstate of L 2 and the expectation value 
of the commutator of x with L 2 over \n£m) vanishes. 
Moreover, L 2 commutes with Hq — u>lL z and therefore 
the sequence of commutators similar the one of Eq. ([63]) 
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can be rearranged so that first one commutes the field op- 
erator with Ho — u>lL z , and next the result is commuted 
with L 2 . However, the expectation value over \n£m) of 
the commutator of any operator with L 2 vanishes, which 
proves our point. Clearly, this result does not depend on 
the orientation of the axis of quantization of the atom 
relative to the external fields. 

Note that even if the time evolution of the expectation 
values of L l and L 2 is only of second order in the external 
field, that does not imply that weak external fields are 
not effective in bringing about changes of the angular 
momentum. For example, in the expansion of Eq. (|63] ) 
the second order terms are multiplied by at least a square 
power of the time; if we consider times comparable to the 
Stark period and the scaling of the matrix elements of the 
position operator with the principal quantum number as 
given in Eq. (28), it is easy to see that the final result is 
not negligible. Moreover, the constraint on the quantum 
invariant % °f Eq. (^) does not say much about the total 
angular momentum of a spherical state, which is not a 
coherent state of the angular momentum, and therefore 
not only one has: 



~y^^(n£m\L t \n£m) 2 ^ (n£m\L 2 \n£m) 



(64) 



but the difference between the two sides of the equation 
can be very large, as one can see by considering a state 
with £ = n — 1 and m = 0. 

Since the results of the previous section show that un- 
der the classical conditions for the external fields inter- 
manifold contributions to the dynamics can be neglected, 
all our considerations apply also to a superposition of 
spherical states with different n's, and therefore our anal- 
ysis sheds some light on the nature of the Rydberg states 
employed in ZEKE spectroscopy. 

In ZEKE, ultrahigh molecular Rydberg states J|-[7|, |||- 
|55| are first excited by a few optical transitions and suc- 
cessively field ionized. This technique is extremely suc- 
cessful because of the ultralong lifetimes of these Ry- 
dberg states, which are explained in terms of extensive 
intrashell mixing of the initial, unstable low-£ states with 
the longer lived high-^ states. For increasing angular mo- 
menta the coupling between the Rydberg electron and 
the molecular core becomes rapidly negligible, so that au- 
toionization and predissociation channels are effectively 
quenched, and the Rydberg state becomes ultralong- 
lived. Therefore it is understood that ZEKE states, i.e., 
the ultralong-living Rydberg states responsible for the 
ZEKE signal, are complicated superpositions of large-n 
spherical eigenstates of the hydrogen atom, which are 
skewed in favor of large angular momentum states. Be- 
cause of the small spacing of high-n Rydberg eigenener- 
gies and of the width of the initial laser pulses, ZEKE 
states initially consist of a superposition of several states 
with different principal quantum numbers |5^,|^). How- 
ever, it is generally assumed that only one angular mo- 
mentum quantum number is allowed in the superposition 



because of the usual selection rules. The population of 
higher-^ states is ascribed solely to the effect of external 
fields. 

In fact, several experimental studies (|-@,[58| have 
shown that the vanishingly small stray fields of the ex- 
perimental set-up and, most importantly, the very weak, 
slowly varying electric fields of the ions present in the 
interaction region populate with great efficacy the high-£ 
Rydberg states which are responsible for the observed ul- 
tralong lifetimes of ZEKE states. On the theoretical side 
recent results 
on the classica 



_6J, |37|, some of which were based 
perturbative approach of Eqs. ^ ]47| - 
|50[ , have explained ^-mixing in terms of the hydrogenic 
model, in which vanishingly small fields are sufficient to 
induce the desired scrambling of the angular momentum 
quantum numbers. The great effectiveness with which 
such extremely weak fields (F < lV/m) populate high- 
l states strongly suggests that the hydrogenic model is 
indeed appropriate to describe angular momentum mix- 
ing in ZEKE states. Moreover, our present findings show 
that the previous classical results [fl7|-|50| are really quan- 
tum mechanical in nature, and can also be extended to 
the case of wave packets. On the other hand, the low- 
l states excited by the laser pulse have a nonnegligible 
quantum defect, which decouples them from the high- 
l, quasi-hydrogenic states; it is then likely that another 
mechanism is at work. More precisely, it is possible that 
the initial optical excitation of the ultrahigh- n states may 
not be strictly limited by the standard selection rules. In- 
stead, by contributions which are of higher order in the 
optical field and yet are nonnegligible because of the ul- 
tralarge dipole moments of Rydberg states -see Eq. ( |2S| ) 
and Eq. (p9[)-, the initial optical pulse may well populate 
a few angular momentum states with relatively large ^'s, 
as one of us has recently shown p7| . Therefore, some 
degree of angular momentum mixing is probably already 
present in the initial Rydberg state, in which case the 
hydrogenic model, in its present extension to superposi- 
tions of states with different quantum numbers, provides 
an accurate description of how weak stray and ionic fields 
bring about the (approximate) randomization in £ and m 
of ZEKE states, which accounts for the observed ultra- 
long lifetimes. 

Finally, the extension of equations Eq. (^) from purely 
classical variables to quantum expectation values lends 
strength to a previous argument of ours concerning slow, 
ion-Rydberg collisions and which until now was based 
solely on purely classical calculations jl7|,fl^]. More pre- 
cisely, we suggested the need for a review of both exper- 
imental and theoretical results for the intrashell transi- 
tions induced in Rydberg alkali atoms by slow collisions 
with ions. In the case of slow ion-Rydberg collisions the 
"magnetic" term of the Hamiltonian arises from the ro- 
tation frequency of the field, and the problem is treated 
in the frame rotating with the field itself. In that frame 
the Hamiltonian is equivalent to the one of a hydrogenic 
electron in weak electric and magnetic fields of constant 
orientation, and time dependent magnitude Jl7|, ff7| ff"9| . 
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The ratio of the two magnitudes, however, remains con- 
stant, and Eq. (§) can be solved exactly @ |9|. It is 
then easy to see that after a full collision, the expectation 
value of L z is over a spherical eigenstate (that is with the 
initial conditions of Eq. (|l|)), is 

46 2 v — 9n 2 

lim (n£m\L z (t)\nim) = —rz — — ; ra (65) 
t-xx Ab^v + 9n z 

where b is the impact parameter of the collision, and v 
is the "reduced" velocity of the incoming ion, i.e. its 
velocity in atomic units multiplied by n, which is the 
principal quantum number of the Rydberg electron in 
the target. In a first approximation, one may insert in 
Eq. ( pq ) an average impact parameter b ~ 50n 2 and a 
reduced velocity v ~ 1, which are consistent with the 
experimental conditions, and the coefficient multiplying 
m becomes > 0.9. Therefore the expectation value of L z 
cannot change much, and if initially the electron is pre- 
pared in a m ^ state, the expectation value of L z will 
not vanish. However, both in the interpretation of the 
experimental data [|62j| . and also in fully quantum the- 
oretical treatments [ |53| , |88|-|9(J the assumption has been 
made of a uniform population of the m substates, which 
corresponds to a zero expectation value of L z . For high 
I states this is a reasonably good approximation, even in 
the case of a non vanishing expectation value of L z . How- 
ever the approximation clearly breaks down for smaller 
values of t, which is precisely the regime for which we 
suggested a critical review of current results. 

V. CONCLUSIONS 

In this paper we have shown that under realistic con- 
ditions, the classical and quantum dynamics of Rydberg 
electrons in weak, slowly varying external fields agree be- 
yond the mandates of Ehrenfest theorem. 

We have shown that for the hydrogen atom in weak, 
slowly varying electric and magnetic fields, to first or- 
der in the applied fields the quantum expectation values 
of the components of the angular momentum and the 
Runge-Lenz vector obey exactly the same equations as 
the time averages (over a Kepler period and along a Ke- 
pler ellipse) of the corresponding classical variables. Our 
proof follows in spirit the approach of classical pertur- 
bation theory, as we fully exploit the properties of the 
zcroth order solutions of the quantum problem, exactly 
as one does in classical mechanics where the time aver- 
aging is done along Kepler ellipses, i.e., the zeroth order 
solutions of the classical problem. Note that this result 
is not an application of Ehrenfest 's theorem, because the 
perturbative approach does not consist of the lineariza- 
tion of the problem in the neighborhood of an equilibrium 
point. Instead, it is an extension, in stronger form, of the 
theorem for the important case of the hydrogen atom in 
weak external fields. 



Most importantly, in our derivation we have not ap- 
plied Pauli's replacement directly in the Hamiltonian 
pl|,|5l],p3[, and therefore we have been able to investi- 
gate the intermanifold contributions to the dynamics. 
In fact, by time averaging the dynamics over a Kepler 
period (which, again, is exactly the same procedure as 
in classical perturbation theory) we have shown that in- 
termanifold terms do not contribute significantly to the 
evolution of the quantum expectation values, as long as 
the strength of the external fields satisfies the same re- 
quirement as in classical mechanics. Interestingly, in the 
semiclassical limit the classical constraint is much weaker 
than the quantum condition for negligible n-mixing, i.e., 
the Inglis- Teller limit. 

This paradox can be resolved by observing that the 
perturbative equations remain accurate only up to times 
comparable to the Stark period, that is, for times which 
in atomic units are ~ 1/ AE, where AE is the energy sep- 
aration between two Stark levels. Therefore, over such 
relatively short times the stationary picture of the Stark 
eigenstates which spread out of a hydrogenic n-manifold 
does not have much physical significance, and that is why 
for these relatively short times the more stringent quan- 
tum condition can be ignored. 

Moreover, by time averaging the intermanifold dynam- 
ics we have also extended the validity of the classical per- 
turbative equations to the case of Rydberg wave packets, 
as long as the spread of the packet over the hydrogenic 
eigenmanifolds is small compared to its average prin- 
cipal quantum number. Note, however, that although 
our analysis shows that the quantum expectation values 
of angular momentum and Runge-Lenz vector evolve in 
time like the classical time averaged variables, it says 
nothing about the localization of the wave packet and the 
quasiclassical dynamics of the packet itself. In fact, our 
time averaging is precisely equivalent to considering a 
spread out version of the wave packet, smeared along its 
orbit. This is the same situation as in classical mechan- 
ics, where one studies the motion of the Kepler ellipse, 
as if it the classical electron had been magically smeared 
along its own trajectory. 

We have also demonstrated that the close quantum- 
classical equivalence can be extended, in the limit of very 
large principal quantum numbers, to the initial condi- 
tions of the equations of motion, provided that the ex- 
pectation values are taken over elliptic states, which are 
states localized along the classical solutions p9|. |(i7|. p8[ . 
Therefore the quantum expectation values of angular mo- 
mentum and Runge-Lenz vector over elliptic states fol- 
low essentially the same trajectories as the time averages 
of the corresponding classical variables. Such complete 
quantum-classical equivalence, however, does not hold for 
the more familiar spherical eigenstates (\n£m)) of the hy- 
drogen atom. 

The realization that the hydrogenic, perturbative 
equations of motion (which account so well for several 
physical phenomena) can also be interpreted as purely 
quantum mechanical equations has led to some insight 
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into the nature the Rydberg states employed in ZEKE 
spectroscopy; it also lends support to our result (previ- 
ously only classical) which indicates that the averaging 
over the L z sublevels (which is used in quantum close- 
coupling calculations to the end of making the problem 
of ion-Rydberg collisions numerically more tractable, and 
also in the interpretation of experimental data) may be 
unjustified. 

Finally, one may wonder if the special equivalence be- 
tween the dynamics of the time averages of classical vari- 
ables and quantum expectation values is a peculiarity of 
the hydrogen atom in weak external fields, or if it can be 
extended to other weakly perturbed integrable systems, 
and the investigation of this problem is in progress in our 
groups. 



Therefore, we consider next the case i = 2, j = 1: 



VPx 



-ly 



Py, L z 



-iyp y L z +i^ y,L z p y + L z yp y } 
-xpy + i^L z yp y ~yp y L z ^ . 



(A3) 



We must then show that the expectation value over \ip n ) 
of the operator within curly brackets vanishes, that is: 



{i/j n \L z yp y - yp y L z \ip n ) = 
Clearly, the state \ip n ) can be written as: 

\ipn) = 22c n (e,m)\n£m) , 



(A4) 



(A5) 



APPENDIX A: PROOF OF THE IDENTITY OF 
EQ. (10) 

In our proof of the special quantum-classical equiva- 
lence of the dynamics of Rydberg electrons in weak ex- 
ternal fields, we have made extensive use of the following 
identity: 



(Al) 



where f t and p 3 are components of the position and mo- 
mentum operator respectively, and where \tp n ) is a state 
confined within a hydrogenic n-manifold. 

In this appendix we prove explicitly the identity of Eq. 
(Al), and we do so for all the pairs of indexes {«, j} to 



stress that our derivation of the equations of motion does 
not depend on the relative orientation between the initial 
axis of quantization of the atom and the direction of the 
applied, external fields in the Hamiltonian of Eq. (||). 
We begin with the simplest case, that is when i = j: 



{^n\XiPi +PtXi\lpn) 

= -i(if) n \x t 



x l: H 



x z ,H 



X t \lpn 



i(tp n \ X^Hq |VO = 



(A2) 



where Hq is the hydrogen atom Hamiltonian and the re- 
sult follows because \4> n ) is an eigenstate of Hq. The same 
approach could be easily extended to all cases. However, 
for the cases in which ?^ja different approach is more 
convenient to the end of studying the intcrmanifold con- 
tributions to the equations of motion, which we do in 
the main text of the paper [see Eq. fl27|)]. Indeed, a 
different proof identifies explicitly the nonclassical terms 
of the Heisenberg equations of motion; these are opera- 
tors which have no counterpart in the classical equations. 
Such terms (see below) yield a null expectation value over 
states which are confined within an n-manifold, and also 
negligible intermanifold contributions to the equations of 
motion (see main text). 



where the C n (£, m)'s are some general coefficients, possi- 
bly complex. By substituting th e ex pansion of Eq. ( A5) 
in the expectation value of Eq. (A.4) one has: 



(ipn\L z yp y - yp y L z \^ n ) 

= ^2^2C n (f, m!)C n {£, m){m'(ril'm'\yp v \n£m) 

&' ,m' l.m / k n \ 

1 (A6) 
— m(n£' m' \yp y \n£m) > . 

where C n (£',m') denot es t he complex conjugate. On the 
other hand, from Eq. (|A2[) it follows that: 



(n£'m'\ypy +p y y\n£m) 

= (n£'m'\2yp y + [p y ,y] \n£m) = , ^ A7 ^ 



and therefore the matrix elements of yp y are: 



(n£'m'\ypy\n£m) = -5 e , e 5r, 



(A8) 



By inserting the matrix elements of Eq. ( A.8) in the dou- 
ble sum of Eq. ( |A6| ) it is easy to see that each term 
within curly brackets vanishes exactly, and therefore the 
identity of Eq. (A4) is proved. 



Next, we consider the case i = 3, j = 1. One has: 

zpx —iz p z , L y 

= izp z L v -ii z,L y p z + LyZp z } 
= -xp z - i \ t v zp z - zp z L v \ , 

and so we must prove that: 

(ip n \L y zp z - zp z L y \ip n ) = 

We use: 



Ly — 



(A9) 



(A10) 



(All) 
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and also: 



L+\n£m) = y/(£-m)(£ + m+ l)\n£m + 1) 
L_\n£m) = y/(i + m){£-m+l)\n£m - 1) , (A12) 



and the expectation value of Eq ( |A10| ) becomes: 

(ipn\L y zp z - zp z L y \ip n ) 



x ly/[t' + m')[£' — m' + l){n£'m' - l\zp z \n£m) 

- sj{£! - m')(£' + m' + l){n£'m' + l\zp z \n£m)^ A13 ^ 

- ^(£~ m)(£ + m + Tj(n£' m'\zp z \n£m + 1) 

+ y/{£ + m){£- m + l)(n£'m'\zp z \n£m - 1)| . 

Clearly, the matrix elements of zp z are also given by Eq. 
(A8), and by inserting that result in Eq. ( |A13 ) it is 
easy to verify that once again the expression within curly 
brackets vanishes. 

Finally, an essentially similar argument proves that the 
identity of Eq. (|Il|) holds also for zp y , which completes 
our proof. 



APPENDIX B: TIME AVERAGING OF THE 
INTERMANIFOLD DYNAMICS 

In this appendix we evaluate explicitly the time aver- 
ages over a Kepler period Tk of the intermanifold con- 
tributions to the equations of motion for the quantum 
expectation values. 

We begin with: 



, T i(E z -E 3 )t\ 



1 



Ik T k 



Tk 



e Ti(E t -E 3 )t dt 



(Bl) 



The integral of Eq. (Bl ) is easily evaluated, and one has: 



,Ti{E t -E 3 )T K _ l 



J_ f K e Ti(E t -E 3 )t dt = ±i l 

Tr Jo \E% ~ EjWk 



However, the energy difference E t — E 3 is: 



(B2) 



r 2 3 r \ j 



(B3) 



where A ZJ = i — j. We then use: 



A, 



j = n 1 + 



(B4) 



where A Jjfi = j — n and n is the principal quantum num- 
ber of the hydrogenic manifold which carries the largest 
weight in the state. The energy difference between two 
manifolds can then be rewritten as: 



o 



A3\ 



where A 3 (i.e., with no indexes) stands for the product 
of any three A's regardless of the indices. The Kepler 
period is: 



T, 



K 



2nn J 



Substituting the results of Eq. ( p5| ) and Eq. (| 
(B2) one obtains: 



t k Jq 
3 



e Ti(E t -E,)t dt 
1 



-A 



O 



A 2 



(B6) 
in Eq. 



(B7) 



Note that to the leading order in A/n the result does not 
depend on the sign of the exponent ; in fact, the leading 
term of the right hand side of Eq. (B7) can be cast in a 
more symmetric form: 



1 1 

Anj ~ 2 _ L ' 1 ~ 2 



1 , 

j,» - n ~ g + 3) 



(B8) 



which concludes the calculation of the first time average. 

Incidentally, by inverting to the leading order the ex- 
pression of Eq. ( p35| ) we obtain a result which we used in 
the main text of this paper: 



l 



Ei — E 1 



(69) 



Next we evaluate: 

'it e -^-E,)t 



K T K J 

Once again, the integral is straightforward: 

— j TK t e-^-^dt = - 1 —— 

Tk Jo [E l — E 3 )Tk 

-i(E % -E 3 )T K _ i 



(BIO) 



{E % - E 3 ) 



(Bll) 



Finally, inserting in Eq. (Bll) the results of Eqs. (B3- 
B6) one obtains: 



TK t i {Et -E 3 )t dt = 

T K Jo A 



1 



6 / 1 



(1 + in A 



hi) 



which concludes our analysis. 
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